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ABSTRACT 

We investigate the local properties of WMAP Cold Spot (CS) localized at (/ = 209°, 6 = -57°) 
in Galactic Coordinate by defining the local statistics: mean temperature, variance, skewness and 
kurtosis. Comparing with the Gaussian random simulations, as well as the other spots in WMAP 
data, we find that the local mean temperatures around CS are quite lower in the scale of < 5° as 
expected. In particular, we also find that the values of local variance and skewness around CS are 
all systematically larger in the scale of > 5°, which implies that WMAP CS is prefer a large-scale 
non-Gaussian structure to a combination of some small structures. The non-Gaussianity of CS is 
totally encoded in the WMAP low multipoles at / < 40. We also find that the cosmic texture can 
excellently explain all the excesses in these statistics. So the local analysis of the WMAP CS supports 
the cosmic texture explanation. 

Subject headings: cosmic microwave background radiation — early universe — methods: data analysis 
— methods: statistical 



1. INTRODUCTION 

Cosmic Microwave Background (CMB) radiation is one 
of the most ancient fossils of the Universe. The obser- 
vations of the NASA Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite on the CMB temperature and 
polarization anisotropics have put tight constraints on 
the cosmological parameters ( Komatsu et ah 2011 ). In 
addition, some anomalies in CMB field have also been 
reported soon after t he release of the WMAP data (see 
(jBennett et al .11201 If ) as a review). Among these, an ex- 
tremely Cold Spot (CS) centered at the position (/ = 
209°, h = —57°) and a characteristic scale about 10° in 
the Galactic Coordinate was detected in the Spherical 
Mexican Hat Wav elet (SMHW) non- Gaussian analysis 
(jVielva et al.ll200l . 

Comparing with the distribution derived from the 
isotropic and Gaussian CMB simulations, due to this 
CS, the SMHW coefficie nts of WMAP data have an 
excess of kurtosis (Cru z et al.l [2005 ). In addition, 
several non- Gaussian statistics, such as the ampli- 
tude and area of the cold spot, the higher criticism 
and so on, have also been applied to identify this 
WMAP CS ([ Cruz et al. 2007a, 2005; Cavon et al>2005"; 
Naselskv et al . 2010; Zhang & Huterer 2010; Vielva 
Since then, various alternative explanations for 
the CS h ave been proposed, inclu ding the po ssible fore- 
grounds (|Cruz et al. 2006; Hansen et al.ll20l"2l ), Sunyaev- 
Zeldovich effect (Cruz et al. 2008), the supervoid in the 
Universe (Inoue & Silk 2006, 2007; Inoue 2012), and the 
cosmic texture ( Cruz et al.l 12007b. 2008). At the same 
time, some analyses have been carried out to distinguish 
these explanations, such as the non-Gaussian tests for 
the different detectors and different f requency channels 
of WMAP satellite (Vielva et al. 200l lCruz et al. 2005), 
the i nvestiga tion of the NV SS sources (Rudnick et al. 
120071 : .Smith fc Huterer II201 0D. the survey around the CS 
with MegaCam. on th e Canada-France-Hawaii Telescope 
(|Granett et aIll2009D. the redshift surve y using VIMOS 
on VLT towards CS (|Bremer et al.ll201Ql ), and the cross- 



correl ation between WMA P and Faraday depth rotation 
map ([Hansen et al.l [20121 ). 

Due to the fact that, nearly all the explanations of CS 
are related to the local characters of the CMB field, the 
studies on the local properties of CS are necessary. In 
this paper, we shall propose a set of novel non-Gaussian 
statistics, i.e. the local mean temperature, variance, 
skewness and kurtosis, to study the local properties of 
the CMB field. By altering the radium of the cap around 
CS, we study the local properties of CS at different scales. 
We find that different from the general properties of the 
foregrounds, the point sources or various local contami- 
nations, in the small scales the local variance, skewness 
and kurtosis values of CS are not significantly large, ex- 
cept for its coldness in temperature. However, after the 
careful comparison with Gaussian simulations, as well as 
the other spots in the WMAP field, we find that when 
R > 5° the local variance and skewness are systemati- 
cally large. This suggests that CS is prefer a large-scale 
non-Gaussian structure to a combination of some small 
structures. In order to confirm it, we repeat the analyses 
adopted by many authors, where the statistics of tem- 
perature and kurtosis in SMHW domain are used. We 
apply these analyses to the WMAP data with different 
^max, and find that nearly all the non-Gaussianities of CS 
are encoded in the low multipoles / < 40. 

It was claimed that the cosmic texture s eems to be 
most promis ing among all the explanations ([Cruz et al.l 
[2007bl [2008[ ) , by investigating the temperature and area 
of CS. In order to check this explanation by our lo- 
cal statistics, we add a similar cosmic texture into the 
Gaussian simulations, and calculate the local quantities 
of CMB fields. We find that the excesses of the local 
statistics of WMAP CS can be excellently explained by 
this non-Gaussian structure. So our local analysis of the 
CS supports the cosmic texture explanation. 

The rest of the paper is organized as follows: In Section 
2, we introduce the WMAP data, which will be used in 
the analysis. In Section 3, we define the local statistics 
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and apply them into WMAP data. In Section 4, the 
non-Gaussianities of CS is studied for the WMAP data 
with different /max, which shows that the non-Gaussian 
signals are all encodes in the low multipoles. Section 5 
summarizes the main results of this paper. 

2. THE WMAP DATA AND SIMULATIONS 

In our analysis, we shall use the WMAP data including 
the QVW7 map, ILC7 map and the NILC5 map. 

2.1. QVW7 

The CMB temperature maps derived from the WMAP 
observations are pixelized in HEALPix format with the 
total number of pixels A/pix = 12A/'g^-^g. In our analysis, 
we use the 7-year WMAP data for Q, V and W fre- 
quency bands with A^side = 512. The linearly co- added 
map (written as "QVW7") is constructed using an in- 
verse weight of the pixel- noise variance cTo/A/obs, where 
(Jo denotes the pixel noise for each differential assembly 
(DA) and A/obs represents the full-sky average of the ef- 
fective number of observations for each pixel. 

2.2. ILC7 and NILC5 

The WMAP instrument is composed of 10 DAs span- 
ning five frequencies from 23 to 94 GHz (jBennett et al.l 
l2QQ3f ). The internal linear combination (ILC) method 
has been used by WMAP team to generate the WMAP 
ILC maps (Hinshaw et al.ll2QQ7l : [Gold et al. 2011). The 
7- year ILC (written as "ILC7") map is a weighted com- 
bination from all five original frequency bands, which are 
smoothed to a common resolution of one de gree. For the 
5-year WMAP data, in (jDelabrouille et al.l l2QQ9) the au- 
thors have made a higher resolution CMB ILC map (writ- 
ten as "NILC5"), an implementation of a constrained 
linear combination of the channels with minimum error 
variance on a frame of spherical called needlet^. In this 
paper, we will also consider both these two maps for the 
analysis. Note that all these WMAP data have the same 
resolution parameter A^side = 512, and the corresponding 
total pixel number A^pix = 3145728. 

In comparison with WMAP observations to give con- 
straint on the statistics, a ACDM cosmology is assumed 
with the cosmological parameters given by the WMAP 
7-year best-fit values (iKomatsu et al.ll2Qlll ): lOOl^b^^ 
2.255, Qch^ = 0.1126, = 0.725, Us = 0.968, r = 0.088 
and A^lko) = 2.430 x 10"^ at ko = 0.002Mpc"^ We 
simulate the CMB maps for each frequency channel by 
considering the WMAP beam resolution and instrument 
noise, and then co-add them with inverse weight of the 
full-sky averaged pixel-noise variance in each frequency 
to get the si mulated QVW7 ma ps. Similar to the pre- 
vious work (jHansen et al.ll2QT2f ). to simulate the ILC7 
maps, we ignore the noises and smooth the simulated 
map with one degree resolution. And for NILC5, we 
co nsider the noise level and beam window function given 
in (jDelabrouille et al]l2009l ). In all the simulations, we 
assumed that the temperature fluctuations and instru- 
ment noise follow the Gaussian distribution, and do not 
consider any effect due to the residual foreground con- 
taminations. 

^ The similar map for 7-y ear WMAP data is recently gotten in 
(|Basak DelabrouilldlMTll V 



QVW7 




Figure 1. The K{R) (left) and K{R) (right) maps for QVW7 
data. In both map, we have adopted R = 2° . 

3. LOCAL STATISTICS OF THE CMB FIELD 

3.1. Local staititics 

In this section, we shall investigate the local properties 
of the CMB field, especially the WMAP Cold Spot, by 
using the local statistics: mean temperature, variance, 
skewness and kurtosis. 

The statistic s of local skewness and kurtosis were firstly 
introduced in (iBernui fc Rebouc as 2009|). For a given 
WMAP data with AT^ide = 512 (QVW7, ILC7 or NILC5), 
we degrade it to the lower resolution N^ide = 256 to re- 
duce the effect of the noises. And then, for this degraded 
CMB map, the constructive process can be formalized 
as follows. Let Vt{Oj^(j)j]R) be a spherical cap with an 
aperture of R degree, centered at {Oj^ We can define 
the functions M (mean temperature), V (standard devi- 
ation), S (skewness) and K (kurtosis) that assign to the 
j*^ cap, centered at {Oj^(j)j) by the following way: 

1 

P i=l 

1 

K,iR) = j^Y.Tt-^, (1) 

^ j i=l 

where A^p is the number of pixels in the j^^ cap, Ti is 
the temperature at i^^ pixel. Obviously, the values 5^ 
and Kj obtained in this way for each cap can be viewed 
as a measure of non-Gaussianity in the direction of the 
center of the cap {6j^(j)j). For a given aperture we can 
scan the celestial sphere with evenly distributed spher- 
ical caps, and build the M(i?)-, V{R)-, S{R)-, K{R)- 
maps. In our analysis, we have chosen the locations of 
centroids of spots to be the pixels in N^ide = 64 reso- 
lution. By choosing different R values, one can study 
the local properties of the C MB field at different scales. 
In (jBernui Reboucas l[2QQ9i . .2010. .20121 the statistics 
Sj{R) and Kj{R) with large R values have been ap- 
plied to study the large-scale globe non-Gaussianity in 
the CMB field. However, in this paper we shall apply 
them to study the CMB local properties. 

It is important to mention that these definitions can- 
not well localize the non-Gaussian sources. For exam- 
ple, in Fig. [T] the kurtosis map K{R) (left panel), we 
can find the clear circular morphology around the point 
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Figure 2. Closewise, M(R) map, V(R) map, S(R) map and K(R) 
map for QVW7, where = 2°. Note that the M{R) and V(R) 
maps have the unit: mK. 



sources. This means that the values of Kj always maxi- 
mize/minimize at the edge of the circles, rather than the 
center of the circles. To overcome this problem and lo- 
calize the non-Gaussian sources, it is better to define the 
following average quantities. 




where X = M, V, K for mean temperature, standard 
deviation, skewness and kurtosis. Xj is the correspond- 
ing local quantities defined above, and A^p is again the 
number of pixels in the cap. For the comparison, we 
plot the corresponding K{R) in the right panel of Fig. 

m 

Now, let us apply the method to the CMB maps. 
Firstly, we consider the QVW7 map. By choosing R = 
2°, we plot X{R) maps in Fig. O The figures clearly 
show that these local statistics are sensitive to the fore- 
ground residuals and various point sources. From M- 
map, one finds that most non-Gaussianity comes from 
the Galactic plane around 6 = 0°. However, from V-, 
S- and iiT-maps, various extra point sources far from the 
Galactic plane are clearly presented. 

Similarly, we also apply the method to the ILC7 and 
NILC5 maps by choosing R = 2°. The results are shown 
in Fig. [3] and Fig. IH We find that these ILC maps are 
much cleaner than QVW7 map in all the four M-, V"-, S-^ 
i^^-maps. Even so, from the S-, iC-maps, we still find 
some no n- Gaussian sources in the Galactic plane. In ad- 
dition, two obvious sources at (/ = 209.5°, b = 20.1°) and 
(/ = 184.9°, b = 5.98°) are clearly localized in ILC7 map, 
which have be identified as the known point sou rces, and 
masked in the KQ75y7 mask ( Gold et al.ll2011f ). 

3.2. The local properties of WMAP Cold Spot 

In this subsection, we shall focus on the local statistics 
for WMAP Cold Spot at (/ = 209°, b = -57°) in Galactic 
Coordinate. First, for a given X{R) map {R = M, V, S 
or K)^ we compare the values centered at CS with those 
of other spots. From Fig. [21 we can find in the maps de- 
rived from QVW7 data, there are too many point sources 
even in the outside of the Galactic plane. So, for a fair 
comparison, in this subsection we shall only consider the 
ILC7 and NILC5 maps. 

To quantify the properties of CS, we define the proba- 




Figure 3. Closewise, M(R) map, V{R) map, S{R) map and K{R) 
map for ILC7, where R = 2° . Note that the M{R) and V{R) maps 
have the unit: mK. 




Figure 4. Closewise, M(R) map, V(R) map, S(R) map and K(R) 
map for NILC5, where R = 2° . Note that the M{R) and V{R) 
maps have the unit: mK. 



bility p- value as follows, 

p = FToh{xi < xo),for M(R), 
p = FToh{xi > xo),for V(R), 
p = Prob(|x,| > |xo|),for S(R),K(R), (3) 

where Xi are the values of X for spots in the region of 
\b\ > 30°, and xq is that for CS. Note that in our analysis, 
the Galactic plane with \b\ < 30° has been excluded to 
reduce the contamination of foreground residuals. So for 
the mean temperature maps, a smaller p- value indicates 
a colder spot of CS, and for the V{R) maps, a smaller 
p- value shows that the variance around CS is larger than 
the other spots. While for S{R) and K{R) maps, smaller 
p- values suggest the larger non-Gaussianities around CS. 
As we can imagine that if CS is simply cold without any 
non-Gaussianity, the p- values for V, S and K should be 
normal, i.e. close to 0.5. On the other hand, if CS is 
a combination of some small-scale non-Gaussian struc- 
tures, as some explanations in (Cruz et al. 2006), the lo- 
cal variance, skewness and kurtosis in small scales should 
be quite large. However, as we will show below, none of 
these is the case of WMAP CS. 

We plot the p- values as functions of R in Fig. [5] for 
both ILC7 (solid lines) for NILC5 (dashed lines). From 
both cases, as anticipated we find that the CS is very 
cold comparing with other spots in WMAP data, i.e. 
p < 10-^ for < 5°. Especially, when = 2°, 3°, 4°, we 
have p = 0, i.e. CS is the coldest spot in the considered 
region of WMAP data. However, as R increases, the p- 
value also increases. In particular, when R > 10°, we 
find p > 0.1, i.e. CS becomes a normal spot. 
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Figure 5. The p- values for ILC7 (solid lines) and NILC5 (dashed 
lines). Black lines are for M, blue ones are for V, red ones are for 
S, and green ones are for K. 



For the V maps, the situation is different: the p- values 
decrease with the increasing of R. When R < 5°, we 
find p > 0.01, i.e. CS is quite normal comparing with 
other spots. But, when R is larger than 10°, p < 10~^ 
is always holden. Similar results are also obtained in the 
S maps. Although it is less significant, we also find that 
local skewness of CS is systematically larger than the 
other spots in the large scale, i.e. R > 10°. Meanwhile, 
Fig. [5] also shows that the local kurtosis is quite normal 
in ah the scales r < R< 15°. 

In order to cross-check these results, we can also com- 
pare WMAP CS with the random Gaussian simulations. 
To do it, we consider 500 random Gaussian samples, and 
calculate the values of X{R) at the spot (/ = 209°, 6 = 
—57°) for each realization. In Figs. [6] and [71 we plot 
mean values and corresponding error bars of the local 
statistics of these realizations and compare them with 
the WMAP results. From these figures, we find the sim- 
ilar results as in Fig. El When < 5°, the WMAP CS 
is quite colder than the simulations. While in the cases 
with > 5°, the local variance and skewness values are 
systematically larger. 

Combining these results, we find that WMAP Cold 
Spot seems to prefer a nontrivial large-scale structure to 
a combination of some small non-Gaussian structures (for 
instance, the point sources or foreground residuals, which 
always follows the non-Gaussianity in the small scales as 
shown in Fig. [2]). This is one of the main conclusions in 
this paper. 

3.3. Comparing with simulation including cosmic 
texture 

In (jCruz et al.ll2007bL [2008'), by studying the temper- 
ature and area of CS, the authors found that the cosmic 
texture, rather than the other explanations, provides an 
excellent interperation for the WMAP Cold SpolQ. In 
this subsection, we shall study if the local anomalies of 
CS found in this work is consistent with the cosmic tex- 
ture explanation. 

^ In the review paper (jVielval [20101 ). one can find some other 
explanations, and the corresponding criticisms, which will not be 
considered in this paper. 
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Figure 6. The values of statistics for ILC7 compare with 500 sim- 
ulations. Red lines denote the WMAP results. Blue ones indicate 
the average values and the corresponding error bars of simulations 
in the case with cosmic texture, and green ones are for the case 
without cosmic texture. 
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Figure 7. The values of statistics for NILC5 compare with 500 
simulations. Red lines denote the WMAP results. Blue ones indi- 
cate the average values and the corresponding error bars of simu- 
lations in the case with cosmic texture, and green ones are for the 
case without cosmic texture. 



The profile for the CMB temperature fluctuation 
caused by a collapsing cosmic texture is given by 



where is the angle from the center, e is the amplitude 
parameter, and i^c is the scale parameter, i^^ = \/3/2'3c- 
By the Bayesian analysis, the texture parameters were 
obtained s = 7^3t|Jx lO^jand i9c 4.9tli^eg at 95% 
confldence (Cruz et al. 2007b). 

In our calculation, we adopt the best-flt texture param- 
eters e = 7.3 X 10~^ and i^c = 4.9°. In above discussion, 
we have compared the WMAP CS with 500 Gaussian 
simulations. Here, in order to taking into account the 
cosmic texture, for each realization we add the texture 
morphology at the same position (/ = 209°, 6 = —57°), 
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Table 1 

For ILC7, the for various statistics 





4 


4 


Xtot 


without texture 35.16 43.27 


31.83 


10.46 


120.73 


with texture 34.01 29.04 


12.76 


7.40 


83.21 


Table 2 








For NILC5, the x for various statistics 
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Xtot 


without texture 33.67 36.23 


43.24 


13.54 


126.68 


with texture 31.17 25.27 


19.82 


11.03 


87.29 



and then repeat the exactly same analyses as in the pre- 
vious subsection. The results are also shown in Figs. [6] 
and [3 (blue dots and error bars). Interestingly enough, 
we find that WMAP CS is exactly consistent with the 
simulations if the cosmic texture is considered. For each 
statistic with every R value, the value of WMAP CS is 
equal to the mean value of simulations in l-a confident 
level. The cosmic texture can excellently account for the 
excesses of M, V and S of WMAP CS. 

In order to quantify it, we calculate the fo^" each 
statistic as follows 



Y{R,))^-'{Yo{R^)-Y{R,)), (5) 



where Y = M,V,S,K. Ri and Rj run through 1° to 
15°. Yo{Ri) are the values of the statistics for WMAP 
CS, and Y{Ri) are those for the simulations. Y{Ri) is 
average value of Y{Ri). S is the covariance matrix of 
the vector Y{Ri). The total value can also be defined as 



Xtot 



For both cases (ILC7 and NILC5) with and without 
cosmic texture, we list the Xy values in Table 1 and 2. 
From both tables, we find that once the cosmic texture 
is considered in simulations, comparing with the pure 
Gaussian samples, every Xy value reduces, especially for 
Y = V and Y = S. So we conclude that the local anal- 
ysis of WMAP CS strongly supports the cosmic texture 
explanation. 

4. COLD SPOT AND WMAP LOW MULTIPOLES 

If WMAP Cold Spot is a large-scale non-Gaussian 
structure as we have found in previous section by the 
local analysis, the non-Gaussianity caused by CS should 
be encoded in the low multipoles rather than the high 
multipoles. In this section, we shall confirm it by study- 
ing the effect of different multipoles on the WMAP non- 
Gaussian sig nals. 

Following dVielya et al.l [2QQ3 : ICruz et al.l [20051 [2QQ6I : 
iZhang Hutererll2QlQD . in this section we study the 
non-Gaussianity of WMAP data by using the wavelet 
transform, which can emphasis or amplify some features 
of the CMB data at a particular scale. The Spherical 
Mexican Hat Wavelets are defined as 

^{0; R) = A(i?)(l + (|)2)2(2 - (lf)e-yy^^\ (6) 
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Figure 8. Values and probabilities of the statistics .Sdisk (upper) 
and Kdisk (lower) for various cases. In all panels, the blue lines are 
for the results of QVW7 map, and red lines are for ILC7 map. The 
circles are for the ILC7 with /max = 20, the crosses are for that 
with /max = 40, and the squares are for that with /max = 100. 



where y = 2 tan(6>/2) is the stereographic projection vari- 
able, and G [0, tt) is the co-latitude. R is the scale, and 
A is the constant for the normalization, which can be 
written as 



AiR) 



27tR' 1 



2 



R^ 
4 



-1/2 



(7) 



The continuous wavelet transform stereographically 
projected over the sphere with respect to ^(^; R) is given 

by 



Tw(r; -R) = y" dn'T{r + f')*(^'; R), 



(8) 



where f = (6>, 0) and f' = iO'^cj)') are the stereographic 
projections to sphere of center of the spot and the dummy 
location, respectively. In our analysis in this section, 
the locations of centroids of spots are chosen to be 
centers of pixels in A^sid e = 32 resolution. Following 
(IZhang fc Huterer |[2Qll . we define the occupancy frac- 
tion as follows to account for the masked parts of the 
sky. 



N^{r',R)= dn'M{r^r')^\0';R), 



(9) 



where M(f) is KQ75y7 mask (|Gold et al.ll2Qll[ ). In order 
to reduce the biases due to masking, we only include the 
results for spot locations r for which N^{r] R) > 0.95. 

In our analysis in this section, we shall consider the 
QVW7 and ILC7 data. We degrade them to a lower res- 
olution with Nside = 128, then apply the KQ75y7 mask. 
For each masked WMAP data, we use the SMHW trans- 
form in Eq. (j8j) to get the corresponding map in wavelet 
domain T^{r;R). To investigate the non-Gaussianity in 
different scales, for each map we consider the cases with 
R = 1°, 2°, 3 °, 4°, 5°, 6° 7°, 8°, 9°, 10° Sim i lar to 
many authors dvTelva et al.l[200l ICruz et al .112 00 5L [20061 : 
'Zhang fc Huterer I l2010 ^.~~we define the statistics 5'disk 
and i^disk as follows to study the non-Gaussianity related 
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to WMAP CS: 

5'disk(-R) - 



^coldest 



(f;i?) 



<Tw(i?) 

1 E':':mih;R) 



spots 



3. 



(10) 
(11) 



Here cr^{R) is the standard deviation of the distribution 
of all spots in a given map, and T^°^^®^*(f, R) is the cold- 
est spot in this distribution. From the definitions, we 
find that 5'disk(^) describes the cold spot significance, 
and i^disk(^) is the kurtosis of spots in a given map. In 
Fig. [8] (left panels), we present the values of Sdisk{R) and 
^disk(^) for different scale parameter R. Both QVW7 
and ILC7 illustrate the same results: the values of both 
^disk(^) and Kdisk{R) maximize at R ^ 5° . The results 
then compare with 2000 randomly generated Gaussian 
simulations, with the exactly same methodology applied. 
So we can get the probabilities of simulations, which have 
the larger Sdisk{R) or Kdisk{R) than those of WMAP 
data. These probabilities for both statistics are also 
shown in Fig. [5] (ri g ht panel). So, similar to other works 
jVielva e t al][200l ICruz et al.l[2005[ iZhang Huto^ 
^010), we find that when = 4° - 6°, WMAP data 
have the deviations from the Gaussian simulations, i.e. 
the corresponding probabilities for the statistics 6'disk(^) 
and/or i^disk(^) are smaller than 1%. 

Now, let us study which multipoles account for the 
non-Gaussianity above. We consider the original ILC7 
map T(f), and expand it via spherical harmonic compo- 
sition: 

T{f) = Y,aimYi„^{f), (12) 

1,171 

where Yim are the spherical harmonics and aim are the 
corresponding coefficients. Then the new map can be 
constructed as follows: 



E 

1=2 m= 



(13) 



It is clear that this new map only includes the low multi- 
poles / G [2, / 

max]- Thus, we can repeat the steps above 
by using the ILC7 map T'{r) instead of T(f). In the anal- 
ysis, we choose three cases with /max = 20, /max = 40 
and /max = 100 to study the effect of different multi- 
poles, and show the results in Fig. [8] with circles, crosses 
and squares, respectively. We find that for the statistic 
5'disk(^), even only the low multipoles / < 20 are con- 
sidered, the values of the statistic and the corresponding 
probabilities are quite close to those gotten in the map 
including all the multipoles. These clearly show that the 
coldness of CS are mainly encoded in these lowest mul- 
topole range, which i s consistent with the conclusion in 
(jNaselskv et al.ll2010[ ). While for the statistic i^disk(i^), 
we find the WMAP data are quite normal for the case 
with /max = 20, comparing with the Gaussian simula- 
tions. However, if /max = 40 is considered, the results 
for both statistics are very close to those in the map 
with all multipoles. So we conclude that WMAP CS 
refiects directly the peculiarities of the low multipoles 
/ < 40, which suggests that CS should be a large-scale 
non-Gaussian structure, rather than a combination of 



some small structures. This is consistent with our con- 
clusion in Section 3. 

5. CONCLUSION 

Since the discovery of non-Gaussian Cold Spot in 
WMAP data, it has attracted much attention, and many 
explanations have been proposed. To distinguish them, 
in this paper we have studied the local properties of 
WMAP CS at different scales by introducing the lo- 
cal statistics including the mean temperature, variance, 
skewness and kurtosis. Comparing with the other spots 
in WMAP data, as well as those in the Gaussian ran- 
dom simulations, we found the excesses of local variance 
and skewness in the large scale at R > 5°, rather than 
in the small scales. At the same time, we found that 
the non-Gaussianity caused by CS is totally encoded in 
the WMAP low multipoles / < 40. These all imply that 
WMAP CS is prefer a large-scale non-Gaussian structure 
to a combination of some small structures. 

It was claimed by many authors that the cosmic tex- 
ture with a characteristic scale about 10°, rather than 
other mechanisms, could provide the excellent explana- 
tion for WMAP CS. By comparing with the random sim- 
ulations including the similar texture structure, we found 
this non- Gaussian structure can excellently explain the 
excesses of the local variance and skewness, in addition 
to the coldness of WMAP CS. So our results in this paper 
strongly support the cosmic texture explanation. 
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